ON THE GEOMETRY OF A DISLOCATED MEDIUM. 



Jeffrey Comer, Ruslan Sharipov 

Abstract. Purely real space versions of the differential equations describing the 
kinematics of a dislocated crystalline medium are considered. The differential geo- 
metric structures associated with them are revealed. 



1. Introduction. 

In [1] a phenomenological approach to the nonlinear theory of plasticity in 
glasses, pitches, and soft polymers was suggested (see also [2]). Currently, we have 
no direct microscopic support for the results of [1] since microscopic mechanisms 
of plasticity in amorphous materials are not yet completely understood, especially 
if one needs an exact quantitative description (see papers [3] and [4] where some 
approaches are developed, but this is by no means the ultimate theory). Having 
no direct way, one should maneuver choosing a detour, a roundabout course to the 
goal. For the theory of plasticity this course goes through the theory of disloca- 
tions (see [5-7]). The matter is that dislocations provide a microscopic mechanism 
explaining the plasticity of crystals. Relying on the integrity of the nature, one can 
expect that the plasticity phenomenon in crystals and in amorphous materials are 
described similarly. 

The paper [5] is a review of the basics. There the nonlinear elastic and plastic 
deformation tensors G and G for a crystalline medium are defined, and the following 
differential equations for them are derived: 



dG 



kq 



dt 



v r v r c kq = -j2 ^kv r c rq - J2 v «« p - 

r— 1 r— 1 r— 1 

3 3 

+ X! ^ + X! ^kr &q, 



(1.1) 



0G[ 
dt 



r=l 

3 



£ v r V r &* = ]T (G[ V r v k - ViV r G r k ) - £ 9 r k G[. (1.2) 



Here v 1 , v 2 , v 3 are the components of the velocity vector v of a point of the medium. 
By V in (1.1) and (1.2) we denote the covariant differentiation. In Cartesian coor- 
dinates the covariant derivative Vj coincides with d/dx l . However, for 
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the sake of generality below we use curvilinear coordinates y 1 , y 2 , y 3 . In this case 
Vj is written through the Christoffel symbols (see [8]) of the standard Euclidean 
metric in the space E (the real space, where the dynamics of any medium occurs): 

ij hi 2 \ 9yi 9yj / ' 

It is remarkable that the same differential equations (1.1) and (1.2) govern the 
evolution of G and G in amorphous media (see [1] and [2]). The difference is that 
the tensorial parameter 9 for amorphous materials is introduced empirically (see 
[1]), while for crystalline materials we have the formula 

3 3 3 

i—l i—1 p—1 

expressing 6 through other tensorial parameters of a medium: T and j (see [6] 
and [7] for more details); T is called the incompatible distorsion tensor and j is 
the tensor of the density of the Burgers vector flow. Unlike 9, both T and j are 
double space tensorial quantities: their upper index i in (1.4) is associated with the 
Burgers space. The tensor field S in (1.4) is expressed through T as the inverse 
matrix: S = T _1 . Its lower index i is associated with the Burgers space. 

The Burgers space B is introduced in [5] as a container for Burgers vectors. It is 
a copy of the real space E filled with the infinite defect-free crystalline grid of that 
material which we have in the real space. The copy of the Euclidean metric g in 
the Burgers space is denoted by g. The Burgers space B is usually equipped with 
Cartesian coordinates . Therefore, we have 



gij = const, g lJ = const . 

Suppose that X is a double space tensor field in E, and assume that i and j are its 
indices associated with the Burgers space: 

.V;:..:;;;;;;; .V;;;.:;;;;;;;// 1 .// 2 .// 1 ;. (1.5) 

Applying S and T to the components (1.5) of the tensor field X, one can convert i 
and j into the real space indices p and q respectively: 

i=l 
3 

X;±;;^;. = J2 x :±:T' f l ( L7 ) 

The elastic deformation tensor G is produced from g according to the recipe (1.7): 



3 3 

Gpq = 9v T p . 



(1.8) 
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The purely real space tensors J and Z are produced according to the recipe (1. 

3 3 



JP q = E § f & Ki = E % V P f r (l-^) 

i=l i=l 

The purely real space tensors R and R are defined similarly (see [7]): 

3 3 

R p q = Y,s?pl i?; 9 = ^sr(v p ^-v ? f;). (1.10) 

i=l i=l 

Here p q are the components of the dual space tensor field p, this tensor field is 
defined as the density of Burgers vector for the dislocations in a crystal (see [5] ) . 
The formulas (1.9) and (1.10) express our intension to write a substantial part of 
the theory in terms of purely real space tensor fields. This goal was declared in 
[7]. The other goal of the present paper is to reveal geometric structures hidden 
underneath the theory of dislocations. 

2. The elastic deformation metric 
and associated connection with torsion. 

Let's consider the elastic deformation tensor G. It is defined by formula (1.8). 
Since det T ^ 0, this formula determines a Riemannian metric in E other than the 
basic Euclidean metric g. Let's call it the elastic deformation metric and denote it 
by G. Then remember the following differential equation derived in [7]: 

3 3 

VpG gA; = 22z£ q G r k +^Zp k G qr . (2.1) 

r—l r—l 

Expressing V p in (2.1) through the partial derivative d/dy p , we get 

3 3 3 3 



dG 



E ~ E ~ E pi *~* rk + E ^p k 

r—l r—l r—l r—l 

Here the Christoffel symbols and T r pk are determined by the formula (1.3). 
Comparing the left and right hand sides of (2.2), we see that it is convenient to 
introduce the other set of Christoffel symbols: 

Ty- = r^- + Z { j. (2.3) 

The Christoffel symbols (2.3) define the other connection in the real space E (dif- 
ferent from the standard Euclidean connection (1.3)) and the other covariant dif- 
ferentiation V. In terms of V the equation (2.1) is written as 

V p G qk = 0. (2.4) 

The equations like (2.4) are known as concordance conditions for metrics and con- 
nections. For instance, g and V are concordant and (1.3) is derived from the 
concordance condition V p g q k = (see [9]). 
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Unlike Iy, the newly introduced Christoffel symbols (2.3) are not symmetric, 
i.e. 7^ r^j. Therefore, they define a nonzero torsion: 

rfik -pk _ -pk In r\ 

± ij ~ L ij ji' K^- ) 

For the torsion tensor with the components (2.5) we use the same symbol T as 
for the distorsion tensor above 1 . However, one should remember that they are 
two different things: the distorsion tensor T is a double space tensor with two 
indices, while the torsion tensor T is a purely real space tensor with three indices. 
Substituting (2.3) into (2.5), for the torsion tensor T we derive: 

K % % (2.6) 

Now let's remember that the tensor fields R and R with the components (1.10) are 
related to each other by the following equality derived in [7]: 

3 3 
s=l r=l 

Here w s y are the components of the volume tensor u) (see [8] and [9] for more 
details). Substituting (2.7) into (2.6), we obtain the equality 

3 3 

^=EE W ^<T< (2.8) 

s— 1 r— 1 

Theorem 2.1. The torsion tensor of the connection (2.3) is determined by the 
density of Burgers vector of a dislocated medium through the equality (2.8). 

Apart from the torsion tensor, any connection possesses another tensorial pa- 
rameter — the curvature tensor. It is given by the formula 

df p df p 3 3 
np o_q m , f m f p _ f m f p (0 Q\ 

T'J — Qyi Q y j 1 im Z^Hqt-jm- K^J 

" m=l ro=l 

The standard Euclidean connection (1.3) is a flat connection, this means that its 
curvature tensor is identically equal to zero: 

dT p dT p 3 3 

f?p ~ n n + pT" r p - \" r m r p = o (i l o) 

dv< dvo J-** oq l im Z^ l *i l o™- v - l z - iu ; 

y y m=l m=l 

Substituting (2.3) into (2.9) and taking into account (2.10), we derive 

%i = V ^0, - V oZ? q + £ Z !rn E Z? q Z P m . (2.11) 

m— 1 m— 1 

1 The matter is that in geometry the symbol «T» is a typical notation for a torsion. 
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In order to calculate ViZ p q and VjZ? q in (2.11) we use the second formula (1.9): 

3 3 3 

V ^J 9 = E V& V.-f ™) = £ S? m V^-f ™ + ]T V/f ™. (2.12) 

rn—l m—1 m—1 

The matrix S is inverse to T. Therefore, we have the equality 

3 3 
n=l fe=l 

Now from (2.12) and (2.13) for ViZ p jq we derive 

3 3 

V l Z p q = Y, S p m V.V.T™ - £ Zf m (2.14) 

m—1 m—1 

and then for VjZf q we write by analogy 

v^r 9 = E & w - - %m K- ( 2 - 15 ) 

m—1 m—1 

Due to (2.14) and (2.15) the above expression (2.11) is transformed to 

Ka = E & (ViV.-f ™ - V,-WC). 

m—1 

Due to the symmetry = 1^ and due to the flatness equality (2.10) the covariant 
derivatives V» and V, are commutative: VjVj = VjVj. Hence, we have 

R p qij = 0. (2.16) 

Theorem 2.2. ^4n?/ dislocated crystalline medium is described by a Riemannian 
metric G (the elastic deformation metric) and by a non- symmetric flat connection 
T being concordant with the metric G. 

The proof is obvious. Indeed, the formula (1.8) provides a metric and (2.3) 
provides a connection. Due to (2.7) this connection is non-symmetric. Due to (2.4) 
this connection is concordant with the metric (1.8) and due to (2.16) it is flat, i. e. 
its curvature tensor is zero. 

The result of the theorem 2.2 is not new. As reported in [10], Kondo, Bilby, 
Bullough, and Smith (see [11] and [12]) in 1950s recognized that dislocations should 
be described in terms of the differential geometry. However, their results are not 
widely known to physicists and engineers. 

3. Reconstructing the distorsion. 



According to the strategy declared in [7], we are going to replace double space 
tensors by purely real space tensorial parameters of a medium and write the com- 
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plete set of the differential equations in terms of these parameters. Suppose for a 
while that this work is done and suppose that the elastic deformation metric G and 
the density of the Burgers vector in its real space form R (see (1.10)) are evaluated 
for some particular medium in some particular case. Then the torsion tensor T 
is also known (see formulas (2.7) and (2.8)). The following theorem says that the 
connection T can be derived from G and R. 

Theorem 3.1. For any Riemannian metric G and for any tensorial field T of 
the type (1,3) there exists a unique connection T concordant with this metric and 
having T as its torsion tensor. 

This theorem is a well-known geometric result. In symmetric case, i. e. if T is 
zero, the connection T is called the standard metric connection or the Levi-Civita 
connection of the metric G (see [9]). 

Proof. The theorem 3.1 is proved by deriving the explicit formula for the compo- 
nents of r. Let's denote by T the symmetric part of T and denote by G the inverse 
matrix for the matrix of the metric tensor G: 

r« = ij n - . Gij = (3.i) 

J 2 

Traditionally, in differential geometry G is called the inverse or the dual metric 
tensor for G and is denoted by the same symbol G. Here we cannot use such 
notations since, apart from G, we have the basic Euclidean metric g, hence, G 1 ^ 
are implicitly determined by the standard index raising procedure: 

3 3 
i=l j=l 

From (2.5) and from the first equality (3.1) we derive 



T ij= T ij + o T *T (3-2) 



Now let's write the concordance condition (2.4) explicitly using (3.2): 

dCr ^ *^ 1 ^ 1 ^ 

~Q~f ~ X! T P<i Grk - F p/c ^r^^E ^PQ Grk + 2 ^ Gq7 " ^'^ 
r—1 r — 1 r— 1 r— 1 

Looking at (3.3), we see that it is convenient to denote 

3 3 

Fkij = ^ I# G rkl T kvj = ^2 T ij G rk- (3.4) 

r—1 r—1 

Since G is a non-degenerate matrix, the transformations (3.4) are invertiblc: 

3 3 
pfc _ \ " p r~irk rpk \ * rp s^rk fo c\ 
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Applying (3.4) to (3.3), wc can rewrite (3.3) in the following form: 

Tfcprj + ^qpk = ~Q~p~ ~ 9 Tkpq — — T qp k- (3-6) 

Performing two cyclic transpositions of indices p — > q — > k — > pin (3.6), we produce 
the other two equalities from the equality (3.6): 

Tgfep + Tpfcq = ^ — - T q kp — — Tpkq. (3.8) 

Now let's add (3.6) and (3.7), then subtract (3.8) from the sum taking into account 
the symmetry of T and the skew-symmetry of T: 

of fdG^,dG kp dG pq \_~ ~ 
kpq ~ { dyP + dyo dy k J pqk qpk ~ 

From this equality, applying (3.5), we derive the following explicit formula for : 

~ k _ x ^ G kr ( dG qr dG rp dG pq ~ ~ \ 

q ~hx 2 \ Q y p °y q d v r qr J' 

And finally, substituting (3.9) into the formula (3.2), we get 

3 



— -s G kr I oGj r oG r i &G 



2 \ dy 1 dyi dy r J 

3 3 S~i rps I (~i rps 

_ ^2 ^2 — — r G kr + - f k 



(3.10) 



ij 



This is the explicit formula for the components of the connection declared in the 
theorem 3.1. Thus, its existence and uniqueness is proved. □ 

A remark. Because of the formulas (2.3) and (2.8) the equality (3.10) is equiv- 
alent to the equality (3.13) from [7]. Thus, the theorem 3.1 yields a geometric 
interpretation for the formula (3.13) derived in the previous paper [7]. 

Now let's return to the equations (1.9). Since S and T are inverse to each other 
(S = T _1 ), the second equation (1.9) is rewritten in the following form: 

3 

y p Tl = Y,Z; q T % r . (3.11) 

r=l 

More explicitly this formula is written as 

Q'J'i 3 3 

^-Ew = E^- ( 3 - 12 ) 

" r=l r=l 
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Applying (2.3) to (3.12), we can bring the equation (3.12) to the following one 

dyP 



(3.13) 



r=l 

And finally, there is the most simple form of the equality (3.13), where V p is used: 

%f l q = 0. (3.14) 

Note that the upper index i does not affect the expansion of the covariant derivatives 
V p and V p in (3.11) and (3.14). This is because the distorsion tensor T is a double 
space tensor and its upper index i is associated with the Burgers space B. 

The equations (3.13) form the so-called complete system of Pfaff equations. The 
section 8 in Chapter V of the thesis [13] and the Appendix A in [14] can be used 
as a brief introduction to the theory of Pfaff systems. All complete Pfaff systems 
are overdetermined systems of partial differential equations. The compatibility 
equations form the basic feature of Pfaff systems. In order to derive them for 
(3.13) one should calculate the second order partial derivatives using (3.13): 

nrpi 3 cips 3 nrfii 3 cvp s 3 3 

01 r _ ST^ UL qr rf^i , -p s 01 s _ UL qr fi ST^ fs pm fi 

dyPdyi ^ dyP i r dyP ^dyv s + ^ ^ « r p s m ' 

y y s=l y s=l y s=l y s =lm=l 

Q-p-^-q = E {^Qyf + E T 7 Km^j t- (3-15) 

Exchanging the indices p and q in (3.15), we derive 

EfS + E^r^fi- (3-16) 



Qyq dyP j-^ \^ dyl 

Subtracting (3.16) from (3.15) and taking into account (2.9), we get 

3 

E ,r n-i ( 3 - 17 ) 

s=l 

The equality (3.17) should be fulfilled for any solution of the Pfaff system (3.13). 
However, in our case we have the stronger result given by the theorem 2.2 and 
the equality (2.16). Regarding the Pfaff equations the equality (2.16) is called the 
compatibility condition of the Pfaff system. 

Let's fix some point Pq within the continuous medium. Without loss of generality 
we can assume that its curvilinear coordinates are equal to zero: y l {Po) = 0. Then 
we choose some constant matrix and denote its components by T's(Q). The equality 

f s i (0,0,0)=fj(0). 

can be understood as the initial value condition for the solution T^y 1 ,y 2 ,y 3 ) of 
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the Pfaff equations (3.13). We prefer to write it in the following form: 

fU =f;(o). (3.18) 

I P=P 

Theorem 3.2. The initial value problem (3.18) for the system of Pfaff equations 
(3.13) has a unique local solution 1 for any predefined matrix T,?(0) if and only if 
the compatibility condition (2.16) is fulfilled. 

The theorem 3.2 is a standard fact of the theory of Pfaff equations. We do not 
give the proof of this theorem here. The idea for its proof can be found in the 
section 8 of the Chapter V in [13]. 

In practice, we need only non-degenerate solutions of the Pfaff equations (3.13), 
i.e. detT ^ 0. Suppose we have two different solutions of the equations (3.13), we 
denote them T[l] and T[2]. Then their initial values at the point Pq are related to 
each other by means of some non-degenerate constant matrix O: 

3 

Tj[2](0)=X)^[l](0)O;. (3.19) 
j=i 

Since (3.13) are linear equations and since the theorem 3.2 provides the uniqueness 
of the solution, the equality (3.19) is fulfilled identically at all points: 

3 

fi[2]=^f|[l]0;. (3.20) 

j'=i 

If both solutions T[l] and T[2] correspond to the same elastic deformation tensor 
G, then from (1.8) and (3.20) we derive 

3 3 
i=l 3=1 

This equality (3.21) means that O is a constant orthogonal matrix with respect to 
Euclidean metric in the Burgers space B. Such a matrix corresponds to a global 
rotation with or without reflection in B. Explaining the concept of the Burgers 
space in [7], we said that it can be understood as an isometric copy of the real 
space E. From this point of view the global rotations and reflections are inessen- 
tial transformations in B. Therefore, up to this inessential uncertainty in defining 
Burgers vectors, now we have the one-to-one correspondence: 

f^G, R. (3.22) 

The double space tensorial field of distorsion T defines the pair of purely real space 
tensorial field: the elastic deformation tensor G and the tensor R, which is the real 
space representation of the Burgers vector density p (see formula (1.10) above). 
For the sake of brevity, from now on, we shall call R the Burgers vector density as 



This means the solution in some neighborhood of the initial point Pq. 
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well. The tensor G is derived from the distorsion T by means of the formula (1.8). 
The formula for tensor R is more complicated: 

3 3 3 3 

R k = E E E E % 9sk <J™ V p f^ . (3.23) 

i— 1 s—1 p—1 q— 1 

The formulas (1.8) and (3.23) correspond to the upper right arrow in (3.22). The 
lower left arrow in (3.22) goes through formula (2.8), through the theorem 3.1 
provided by the formula (3.10), and through solving the system of Pfaff equations 
(3.13). As a result we recover the incompatible distorsion tensor T from G and R. 

A remark. Being produced from the single field T, the tensor field G and R 
are not absolutely independent. They are related to each other through the zero- 
curvature condition (2.16). 

4. Time evolution and consistence 
of the kinematic equations in whole. 

The time evolution of the elastic deformation tensor G is given by the equation 
(1.1). The time evolution for the tensor R is given by the equation (4.11) from [7]: 

dR k 3 3 33 

-gf - E J P k R l - E v m v k i^-EE vP zt P k + 

p—1 rn— 1 rn— 1 p—1 

3 3 3 3 3 3 3 \ ' ) 



+ E E E 9*» " mrs ^ J k s + E E E E ^ " mra z% j p s = o. 

m—l r—1 s—1 m— 1 r—1 s—1 p—1 

Now let's remember the second relationship (1.9) and apply it to (3.23). Then 

3 3 3 

^=EEE^ wmrs ^ s - ( 4 -2) 



m—l r—1 s—1 



Substituting (4.2) into the second term of (4.1) and rearranging the terms, we get 



dR f - E ( v '» vfe + E vP ) R * 



at 

m—l \ p—1 



(4.3) 



+ E E E 9,m US™* UrJ k s + J] Zr P ^ ~ 4 ) = 0. 

m—l r—1 s—1 \ p—1 / 

Then we remember the relationship (2.3) and replace V by V in (4.3): 

dR k 3 3 3 3 

]T V m v k R™ + E E E B<m " mrs V r J k = 0. (4.4) 



dt 

m—l m—l r—1 s—1 



In a similar way let's replace V by V in the evolution equation (1.1) for the elastic 
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deformation tensor G. Applying (2.3) to (1.1) and using (2.4), we derive 
^ + £ V fc «' G rq + £ £ «<" - Zf r ) G P9 + £ V &r + 

r— 1 r— 1 p— 1 r— 1 

+ E E yr (%* - fy) d *> = E ^ + E 

r— 1 p— 1 r— 1 r— 1 

Now let's compare (1.9) with (1.4). As a result we write (1.4) as 

3 



(4.5) 



P =i 

Applying (4.6) to (4.5), we transform (4.5) to the following equation: 
^ + J2 Vkv r G rq + J2 V q v r G kr + 

r—1 r—1 
3 3 

+ E J k Gr 1 + E Jq = 0. 

r—1 r—1 



(4.7) 



The differential equations (4.4) and (4.7) describe the time evolution of the tensor 
fields G and R. The next step is to show that this time evolution is compatible 
with the zero-curvature condition (2.16). The following commutation relationships 
are derived from the definition of covariant derivatives by direct calculations: 

3 opfe 

[du v p] x k = Y,^r xs > ( 4 - 8 ) 

S — l 

3 gf, s 

[ dt ,V p ]X k = -J2-^Xs- (4-9) 

s—l 

Here dt = d/dt, while X k and Xj. stand for the components of arbitrary vectorial 
and covectorial fields respectively. In order to calculate the time derivatives of the 
connection components in (4.8) and (4.9) we use the equality (2.4) rewritten as 

dC 1 

" E f ^ ~ E n k G qs = 0. (4.10) 

y s=i s=i 

For the sake of brevity in the further calculations we denote 

dr^s _ p k dGrs _ a dT^s _ ^ 

dt rs ' dt ™ rs ' dt 

where T k s are the components of the torsion tensor T related to R through (2.8). 
Now differentiating (4.10) with respect to t, we find 

3 3 

Y,hAk+T,hkG gs =%(G gk ). (4.ii) 



12 



JEFFREY COMER, RUSLAN SHARIPOV 



The equality (4.11) is similar to (3.3). Therefore, we shall treat it similarly. Let's 
apply (3.1) and (3.2) to (4.11). As a result we get 



E r p i sk + E r p k I s — 

s=l s=l 
X 1 3 X 1 3 X 

%{G q k) - 2^2^pq^ sk ~ 2 E^p**"*' 8- 

8=1 8=1 



(4.12) 



Now we introduce the following notations similar to (3.4): 
3 . 3 

f k ij = E r y Grk , T k ij = ^ Tij G rk ■ (4-13) 

r— 1 r— 1 

Applying (4.13) to (4.12), we strengthen the resemblance of (4.12) and (3.3): 

f kpq + F qp k = V p (G g fe) — — Tkpq ~ — tqph- (4-14) 

Note that (4.14) looks pretty like the equality (3.6). Therefore, we can use the 
same arguments as in proving the theorem 3.1 and derive the following formula: 



h = E ?(Vi(4>) + %(G n ) - v r (6 tJ ) 



2 

r=l 



(4.15) 

^ jr T" ^js * ir gfc r , * j* 



- 1 ,> T Ujj i ,;„ ~. 1 ~, 



r— 1 s— 1 



The covariant derivatives Vi(Gj r ), Vj(Gh), and V r (Gy) arc given by the formula 
(4.7). Applying this formula, we get 



Vi(Gjr) + Vj(Gri) + V r (G«) = - £ (V ? V, V m + 

m— 1 

3 

■ V,- V, « m + V, J™ + Vj Ji m ) G mr + ([Vr, V]u ro + V-JT - (4.16) 



VJ™) G mj + ]T ([V P , %]v m + V r J™ - Vj J, m ) G 



m— 1 



As for the time derivatives T s r , T| r , and Ty, we use the formula (2.8) differentiating 
it with respect to t. As a result, applying (4.4), we obtain 



3 3 3 > 

k 



f % = E E ^ -9 sr E - E E E »™ w " Q/3 ^ J / 

r— 1 s— 1 \ ?i— 1 n.— l a— 1 /3— 1 
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The above huge formula can be simplified to the following one: 

3 

^-V^-V^ + E^y^- (4.17) 

n=l 

From (4.17) we easily derive the following two equalities: 

3.3/ 3 \ 

£ G ia T s jr - £ G im V r J™ 1 - VjJ r m + £ V n v m Z£ , (4.18) 

s— 1 rn—1 \ n—1 / 

3 3 / 3 \ 

£ G js Z? r = £ G jm V r Jr - V, J™ + £ V„ W " 1 Z™ . (4.19) 

s— 1 rn—1 \ n—1 / 

Now we apply (4.16), (4.17), (4.18), and (4.19) to (4.15). As a result we obtain 
^ = 2 J V,J/ + - £ V„^ T« + 

m— 1 r— 1 \ ri — l / 

m— 1 r— 1 \ n—1 / 

The commutators [V r , Vj] and [V r , Vj] in the above formula (4.20) can be calculated 
on the base of the well-known differential-geometric formula 

3 3 

[V i; Vj]X k = £ R% X n - £ f » V n X k . (4.21) 

n—1 n—1 

Here X k and X n stand for the components of an arbitrary vector field, while i?„ _ 
are the components of the curvature tensor given by the formula (2.9). Choosing 
X = v in (4.21) and substituting (4.21) into (4.20), we derive 

~ "T 1 v (4.22) 

3 3 3 I n kr r>m „,n fi nkr f>rn ,.n n . \ 

_|_ / Jri -nri v ^mj u J Hirj u j 

m— 1 r— 1 n—1 \ / 

Now we are ready to calculate the time derivative of the curvature tensor R due 
to the time evolution of tensor fields G and R given by the equations (4.4) and 
(4.7). Differentiating (2.9) with respect to t, we get 

dR k ■ ^ 

= Vif* - %ff q + £ Z S*f* (4.23) 

rn—1 
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Before substituting (4.22) into (4.23) we transform it using the identity (4.21) again: 

3 3 

% % v k = WiWj v k -Y, Rnij v n +J2 f % V n v k - 

n—1 n—1 

Due to this identity, the equality (4.22) now looks like 

1 3 



2 

n=l 



3 3 3 / ^,fc r A m „ A (~,kr nm n n 

/ Or rt nri V G mj - 

m— 1 r— 1 n—1 V 



Continuing our calculations, below we shall omit the terms containing the compo- 
nents of the curvature tensor and their spatial derivatives denoting them by dots: 

t% = v k - V, J k + . . . . (4.24) 

Then, substituting (4.24) into (4.23), we derive: 
dR k 



at 



= -ViVjVg v K - ViV,- J q K + Vj-ViV, v 



^AJ q k - E % v™% yk - E ^ m J < 



(4.25) 



fc 

m— 1 



In order to simplify this formula let's write (4.25) as follows: 
— ^ = -[Vi, V,-] (%^ + J £ fc 



3 

-E% v m (y q v k + j q k 



(4.26) 



The further transformation of (4.26) is based on the other well-known formula of 
the differential geometry, which is similar to (4.21): 

3 3 3 

[% %]X k = ]T R k nij X n q - £ R n qij K V n X k . (4.27) 

n — 1 n—1 n—1 

Here X k stand for the components of an arbitrary tensorial field of the type (1,1). 
Substituting X k = V q v k + J q k into (4.27), we see that (4.26) can be written as 

dR k 3 3 

= - E k ni3 (% + J,") + E (^n V " + J n) + ■ • ■ ■ ( 4 ' 2 8) 

n—1 n—1 

Looking at (4.28), we see that two terms explicitly written in the right hand side 
of this formula contain the components of the curvature tensor R. According to 
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the above our convention, they should also be denoted by dots. This means that 
all terms in the right hand side of (4.28) would vanish provided the zero-curvature 
condition (2.16) is fulfilled. 

Theorem 4.1. The time evolution of the tensor fields G and R given by the equa- 
tions (4.4) and (4.7) preserves the zero-curvature condition (2.16), i. e. the curvature 
tensor R is permanently equal to zero if R = at some initial instant of time. 

Passing from T back to the standard Euclidean connection T in the real space 
E, and hence, from V back to V, we can formulate the theorem 4.1 as follows. 

Theorem 4.2. The time evolution of the tensor fields G and R given by the equa- 
tions (1.1) and (4.1) preserves the zero-curvature condition (2.16), i. e. the curvature 
tensor R is permanently equal to zero if R = at some initial instant of time. 

5. Zero-divergency condition. 

The zero-curvature condition (2.16) is not the only condition the tensor fields 
G and R should obey. Another condition was derived in [7] from the equality 
div p = 0. Therefore, we call it the zero- divergency condition: 

jzjZa pq v P R k q + E E E s pq K z k vm = o (5.i) 

p— 1 q— 1 m—1 p— 1 q— 1 

(see (4.8) in [7]). The following theorem shows that (5.1) is not an independent 
condition for G and R. 

Theorem 5.1. The zero-divergency condition (5.1) can be derived from the zero- 
curvature condition (2.16). 

Formulating the theorem 5.1, we assume that the tensor fields G and R are given 
and that the following conditions are fulfilled: 

(1) R define the torsion tensor T according to the equality (2.8); 

(2) G and R define the non-symmetric connection T with the torsion T according 
to the theorem 3.1; 

(3) the curvature tensor R of the connection T is equal to zero; 

(4) the tensor field Z is determined as the difference of the connection T and the 
standard symmetric Euclidean connection T according to the formula (2.3). 

From the above four conditions (l)-(4) one easily derives the equality 

3 3 3 
n— 1 r—1 s—1 

and the formula (2.11) for the components of the curvature tensor R. The zero- 
curvature condition R = then is equivalent to the symmetry condition 



U k = U k 

^ spr ^ s i 



(5.3) 
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where the quantities U^ pr are given by the formula 

3 

Ugpr = ^P^rs + E Z™s Zpm- (5-4) 

m— 1 

Due to the symmetry condition (5.3) the following sum is obviously equal to zero: 

EEI> prs[ C = o. (5 - 5) 

p—l r—l s—1 

Indeed, U^ pr is symmetric, while u> prs is skew-symmetric with respect to p and r. 
The next transformation of the identity (5.5) is also quite obvious: 

E E E E E 9 P q a qn w - rs < P = o. (5.6) 

g— 1 n— 1 p—l r—l s—1 

Now it is sufficient to substitute (5.4) into (5.6). As a result we obtain 

3 3 3 3 3 



EE^EEE^-v p i 

p—l q—1 71 — 1 r—l S—1 



k 

r s 



33 3/333 \ 

EE^ 9 E EEE^""t z 

p—l q—l rn — 1 \ n— 1 r—l s—1 / 



(5.7) 



i 0. 



Applying (5.2) to (5.7), we find that the equality (5.7) is equivalent to the zero- 
divergency condition (5.1). Thus, the theorem 5.1 is proved. 



6. Conclusions. 

As a main result we can formulate the following statement: the elastic deforma- 
tion tensor G and the tensor of Burgers vector density R are two basic tensor fields 
describing completely the deformation state of a crystal. The time evolution of G 
and R is determined by the differential equations (4.4) and (4.7) 1 . However, (4.4) 
and (4.7) do not form the closed system of differential equations — they are only 
the kinematic equations. They should be completed with the dynamic equations 
relating J (the density of the Burgers vector flow) to G and R. Qualitatively, the 
process of a crystal deformation is expressed by the following diagram: 



R, G 


i 


a 


2 


J 


> » 


J » 





J 



The elastic deformation G produces the stress <x, this phenomenon is expressed by 
the arrow 1 (Hooke's law or its nonlinear generalization). The stress ex causes the 
dislocations to move producing their flow J, see the arrow 2 on the above diagram. 



1 Or by the equations (4.1) and (1.1), which are equivalent to (4.4) and (4.7). 
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And finally, the moving dislocations rearrange the interatomic bonds causing G and 
the stress <x to relax. This phenomenon is expressed by the arrow 3 and described 
by the equations (4.4) and (4.7). The plastic deformation tensor G is not presented 
on the diagram. However, it can be implicitly present in the arrows 1 and 2. The 
detailed quantitative description of the phenomena associated with these arrows is 
the subject of the separate paper. 
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